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Abstract 



We examine the time-reversal-violating nuclear "Schiff moment" that in- 



duces electric dipole moments in atoms. After presenting a self-contained 
derivation of the form of the Schiff operator, we show that the distribution of 
Schiff strength, an important ingredient in the ground-state Schiff moment, is 
very different from the electric-dipole-strength distribution, with the Schiff mo- 
OS ■ ment receiving no strength from the giant dipole resonance in the Goldhaber- 

Teller model. We then present shell-model calculations in light nuclei that 
confirm the negligible role of the dipole resonance and show the Schiff strength 
to be strongly correlated with low- lying octupole strength. Next, we turn to 
q heavy nuclei, examining recent arguments for the strong enhancement of Schiff 

moments in octupole-deformed nuclei over that of 199 Hg, for example. We 
concur that there is a significant enhancement while pointing to effects ne- 
glected in previous work (both in the octupole-deformed nuclides and 199 Hg) 
that may reduce it somewhat, and emphasizing the need for microscopic calcu- 
lations to resolve the issue. Finally, we show that static octupole deformation 
is not essential for the development of collective Schiff moments; nuclei with 
strong octupole vibrations have them as well, and some could be exploited by 
experiment. 



1 Introduction 

Observation of an atomic electric dipole moment would signal the violation of time- 
reversal (T) symmetry [Q, which kaon decay tells us is present at some level[||]. So far all 
measurements, whether on elementary particles or atoms and despite rather high sensi- 
tivity, have been statistically zero, but experiments continue to improve [[|. The level at 
which dipole moments are finally seen will help decide among a number of candidates for 
the fundamental source of T violation. 

Several theorists have proposed that the light actinides would be the best elements 
in which to detect a small dipole moment[Q, [B|, ^. Most recently, the authors of refs. 
P, pi have argued that the existence of octupole (pear-shaped) deformation in the nuclei 
of these atoms enhances the sensitivity of atomic dipole moments to nuclear parity (P) 
and T violation by factors of 100 to 1000 (typically about 400 in the later reference) over 
the sensitivity in the atom with the best current experimental limit, 199 Hg. This level 
of enhancement is due in large part to the existence of close-lying parity doublets and 
favorable atomic structure in the light actinides, but also to fact that it is not the dipole 
moment of the nucleus that induces a dipole moment in the surrounding electrons, but 
rather the "Schiff moment" , a quantity that reflects the mean-square radius of the nuclear 
dipole distribution. Asymmetric nuclei have large intrinsic Schiff moments even though 
their intrinsic dipole moments are very small, in the same way that a neutral particle can 
have a finite charge radius. 

The arguments of refs. P, warrant careful investigation. In this paper, we give a ped- 
agogical derivation of the Schiff operator, explore its action on nuclear ground states, and 
address the role of octupole correlations in generating ground-state Schiff moments. The 
discussion is organized as follows: section I contains a derivation of the nucleus-electron 
interaction responsible for atomic-dipole moments and introduces the Schiff operator. The 
section concludes with an evaluation of the nuclear Schiff moment under the assumption 
that the dominant source of time-reversal invariance is a nucleon electric dipole moment. 
Section II reveals important differences between the dipole and Schiff operators, showing 
that in the Goldhaber- Teller model no Schiff strength is produced by the giant dipole 
resonance. In section III, we look at Schiff moments in light nuclei, particularly 19 F, 
confirming the near absence of Schiff strength in the giant resonance and pointing out 
a strong component of strength correlated with low-lying octupole excitations. The first 
excited state provides the largest contribution to the Schiff moment. Section IV takes up 
octupole correlations in heavy nuclei, focusing on octupole deformation. We find no flaw 
in the argument that moments in such nuclei are collective and enhanced, but point to 
physics that may make the enhancement less dramatic than claimed in refs. @, ^ (more 
detailed microscopic calculations of both the octupole-deformed nuclei and those that are 
currently used in experiments should resolve the uncertainty) . In Section V we argue that 
the collective Schiff moments do not depend on the delicate and sometimes unanswerable 
question of whether a nucleus is octupole deformed. Low-lying octupole vibrations gener- 
ate them in the same way as static octupole deformation, increasing the number of atoms 
in which one can expect large effects. Section VI summarizes our findings. 



2 Schiff moments 

We begin by deriving the nucleus-electron interaction responsible for generating atomic 
dipole moments. Though the result is well known || [5], |j, a complete derivation has never 
appeared in one place, and our derivation differs in its details from the others. 

In 1939 Feynman|| developed a quantum theory of molecular forces as part of what is 
now called the Hellmann-Feynman theorem, or sometimes the parameter theorem. This 
charming relative of the virial theorem|lC] allows insight into how forces in a complex 



quantum system are balanced against one another, just as they are in a classical system. 
In a molecule, for example, internal Coulomb forces between electrons and nuclei counter 
each other so that there is no net force on the molecule (or it would move). In the system 
we consider here, a neutral atom in a uniform electric field, there is no net charge so there 
is no net force on the system. To achieve this the electrons rearrange themselves [11] so 
that there is no net electric field at the nucleus (or it would move). 

This shielding effect has dramatic and unfortunate implications for experiments that 
would probe the atom-nucleus system. Because nuclei are of finite extent, however, the 
shielding of the electrons varies over the nuclear volume, and this makes probing the 



nuclear interior possible, as shown originally by Schiff [12]. We present next a variant of 

Schiff's derivation that uses modern effective- field-theory techniques fi~3|| to produce the 

simpler approximate result that has been obtained more recently ffl, Rl H. 

We assume a neutral (nonrelativistic, for simplicity) atom containing an extremely 

heavy nucleus with non- vanishing spin sitting in a uniform electric field, Eq. The atom 

contains Z electrons, each with charge e, while the nuclear charge is Ze p , where e p = — e is 

the proton charge, and in our units the fine-structure constant is given by a = ei/Ai:. The 

nucleus has both an electric monopole distribution and a tiny electric dipole distribution 

leading to an electric dipole moment d = e p do; other moments can be easily added. In 

electric-field gauge the atomic-plus-interacting-nucleus Hamiltonian can be written in the 

form 

z 

#atom = ^2 [ Ki + Vi + e< K^) - eE ■ n - e p E ■ d , (1) 

i=i 

where -JQ is the kinetic energy of the i— electron (pf /2m e in the nonrelativistic approx- 
imation), r*j and pi are that electron's coordinate and momentum relative to the nuclear 
center-of-mass (CM), m e is the electron mass, 

^=«Ep^i ( 2 ) 

is the electron-electron Coulomb interaction, and 



i „ f d 3 x p(x) 

4tt 



P _ , ^j^_ (3) 



is the electrostatic potential due to the complete nuclear charge distribution p(x) (with 
dimensions ^~ 3 , and normalized to Z). 



To explore the physics of eq. (|l]) we employ a trick to remove the last term in that equa- 
tion. Because |do| = ^o is tiny even on nuclear scales, it is sufficient to manipulate eq. (1) 
to first order in that quantity, ignoring all higher-order terms. For example, performing a 
first-order unitary transformation on H a tom produces i/ atom — -ffatom + i [U, -ffatom] , where 
U ~ do m our case is t ne Hermitian operator 

CT = |-£tf. (4) 

Performing the commutator generates two terms, one from the last bracketed term in 
eq. (1) that cancels the interaction of the nucleus and the external field (U was constructed 
to do this), and a second that takes its place: 

Z , s 

H a tom = Y,( K i + V i- e ^o-U + eHn)-^d -V t ^(f i )\ . (5) 

The difference, A, between eq. (5) and eq. (1), which cannot lead to an energy shift to 
first order in do, is given by 




A = e p d ■ [E - - \ V, <P(fi) = . (6) 



This type of relationship, generated by obvious equalities such as ([U, -Hatom]) = 0, is 
often called a hypervirial theorem[^| and can be derived from the Hellmann-Feynman 
theorem^. If one neglects the finite extent of the nucleus and replaces p{x) by Z5 3 (x), 
eq. (6) can be rearranged into the form 

Apoint = e p d ■ (E + E e ) = , (7) 

where E e is the electric field at the nucleus caused by the electrons. This simple result 
states that exact screening holds in the point-nucleus approximation, and highlights how a 
nonzero nuclear volume leads to a small but significant breakdown of screening. Screening 
is now directly incorporated into H at om- 

We can now take advantage of the two very different scales -- atomic and nuclear 

- in the Hamiltonian. All of the nuclear physics is contained in p(x) and reflected in 

4>{r); since (Rn/Ra) ~ (lfm/1 A) ~ 1CP 5 , only a few moments of p(x) will have practical 

importance. For this reason we apply a derivative expansion of the type used in effective- 



field theories |13| to p(x). We assume that the monopole part of p(x) can be expanded in a 
series of the form: a <5 3 (x) + b V 2 J 3 (x) + • • •■ An analogous expression holds for the dipole 
part. For this to make sense the coefficients a,b, . . . together with the derivatives must 
reflect increasing powers of Rn/Ra- Moreover, we must preserve conventional definitions, 
such as 

d 3 x p(x) = Z , (8) 



d 3 xx 2 p(x) = Z(r 2 ) ch , (9) 



d x x p(x) = do 
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(10) 
(11) 



where the vector quantity On [15], which is the second moment of the dipole distribution, 
bears a similar relationship to do as Z(r 2 ) c h does to Z. Thus, we posit 



p(x) 



Z5 3 (x) + Z^V 2 5 3 (x) 
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(12) 



as a sum of monopole and dipole parts. Because derivatives with respect to x in (j) (viz., 
from p) can be transformed into derivatives with respect to r*j through integration by 
parts, this is indeed an expansion in Rn/Ra- 

We next separate -Hatom into a part independent of dipole moments do and Oq and 
another time-reversal- and parity-violating part proportional to these moments: 



•" atom 

Expanding p(x) as in eq. (|45j) , we get 
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\x — f\ Z J \x — r\ 

Writing out the explicit expansions for p mon and paip leads to the general result for Ah 
expressed in terms of the Schiff moment B, S: 



Ah(f) 



4TrS-V5 3 (r) + - 



s = 10 



-> 5 - 
Oo ~ ^d (r 2 ) ch 



(17) 



Thus, the coupling of the nuclear dipole distribution to the atomic electrons is through the 
Schiff moment^. The result ( |T7| ) depends in leading order on terms of order R 3 ^, because 
eq. (0) mandates the cancellation of terms of order Rjy. 

One can make contact with Schiff's paper[12] by defining quantities pc(x) and Pm(x) 
such that p mon (x)/Z = pc{x) and Pd\ p {x) = —do ■ V pm(x). Equation (10) of ref. [[L^] then 
follows from eq. ( |l~6| ) above: 



Ah 



Schiff/ 



do -f 



d 3 x(p M {x)- p c {x))6{r -x) . 



(18) 



1 The factor of A-k in the first of eqs. IX7| is often H, W] incorporated in the definition of S. 



This result is exact but not particularly useful. Expanding pc and pm in the form of eq. 
(|l~2l) leads to the approximate result 

Ah SchiS (r) ~ yd • V5 3 (r) (<r 2 ) M - (r 2 ) c ) , (19) 

from which we deduce his form] 2 ] of Oq: 

dl chiS = 5 -d (r 2 ) M . (20) 

Finally, the hypervirial (or Hellmann-Feynman) quantity A in eq. (6) can be written 
using eq. (12) in the form 

A = e p d -(E + E e + V 2 E e ^-^ + ---) = , (21) 



where E e is the electrons' electric field at the nuclear CM, and the last term arises from 
averaging that electric field over the nuclear volume. Because it is the averaged field that 
cancels the external field at the nuclear CM, we see explicitly how nuclear finite size affects 
screening. Moreover, because that last term is equivalent to 

y^E^'^^) . ( 22 ) 

we can see that it produces the do term in the Schiff moment in eq. (|i~7|) (via the second 
term in eq. fll6|))- The Oq term in eq. (O) arises from the first term in eq. (|i~6|). 

Having formulated expressions for nuclear Schiff moments, we want to use them to- 
gether with assumptions about the dominant source of P and T violation to evaluate Schiff 
moments in real nuclei. In the rest of this paper we will assume that a P,T- violating com- 
ponent of the nucleon-nucleon interaction causes a Schiff moment in the distribution of 
protons, but we conclude this section by briefly describing another possibility: that dipole 
moments of individual nucleons are responsible for the nuclear Schiff moment^. 

We introduce proton and neutron (isotopic) projection operators, pi and hi, for the i— 
nucleon. The dipole moment must point along the spin, <?j, of the nucleon, leading to the 
impulse-approximation result (7| 

Pdi P (r) = (^PidpBi-Vi p p PT (fi -r)+hid n ai- V* p% T {n - r) \ , (23) 

where p PT and pp T are the proton and neutron electric dipole densities (normalized to 1) 
associated with d p and d n . This yields 

d "o = ( J2 (Pi d P + *i d n) °i) = 4 + do , (24) 



2 This result is a little misleading because it seems to imply that the magnitude of Oo should be do 
times a typical nuclear size; this need not be the case. 

3 Meson-exchange currents can also generate P,T-violating nuclear moments directly. 



which depends only on the ground-state expectation value of nucleon spin and isospin 
operators. The Schiff moment can be obtained by evaluating Oq with eq. (J23[), producing 

^ = f {(r 2 ) P PT + (rVp t) + f {(r 2 )PT + (tTpt) , (25) 

where (r 2 ) p PT and (r 2 } PT are the mean-square radii of the densities pp T and pp, T and 

r£Pid P ^rf\ = 4(r 2 ) z PT ; fe^e^rA = cT^(r 2 )£ T . (26) 

We expect (r 2 )p T to be comparable to nuclear sizes, and thus much larger than (r 2 ) 1 ^, 
which should be comparable to nucleon sizes. 

3 Distribution of Schiff strength: its significance and the 
role of the giant dipole resonance 

For our considerations, as we have said, the most important mechanism for a nuclear 
ground-state electric dipole or Schiff moment is the action of a pseudoscalar T-violating 
nucleon- nucleon potential, Vpp- This induces a Schiff moment given by second-order 
perturbation theory: 

SS ft)=E WI !'"f" W + c.c, (27) 

where the state |^o) has J z = J ^ and S z is the z-component of the Schiff-moment 
operator S. The distribution of Schiff strength to the excited states |\Pj) is therefore 
a crucial ingredient in the ground-state moment. At first sight, one might think that 
this distribution should resemble that of the electric dipole operator fr z , but the two are 
remarkably different. Most of the electric dipole strength is in a broad resonance at 10 
to 20 MeV of excitation energy. Almost none of the Schiff strength, however, goes to the 
states in the giant resonance. In fact, in the simple but venerable Goldhaber- Teller (GT) 
model, in which the giant dipole resonance corresponds to the \Tiui oscillation of all the 
protons with respect to all the neutrons, the Schiff strength to the resonance is identically 
zero. 

To see this, we start by assuming [jig] rigid distributions for the Z protons and for the N 
neutrons, which oscillate harmonically about their CM with frequency to. The separation 
of the two rigid spherical (for simplicity) distributions is denoted q and the distance of the 
proton CM from the overall CM is Nq/A. The charge distribution is then given by 

PGT(x)=p (\x-Nq/A\) , (28) 

where po is the "bare" nuclear- ground-state charge distribution normalized to Z protons. 
Expanding this equation, one finds components of /Sgt proportional to \q\, \q\ 2 , \q\ 3 , •••, 



or equivalently (because q oscillates harmonically), Vfuu, and 2huj, 1 and 3Huj, etc. In 
particular one finds that the g 2 -term renormalizes the ground-state charge distribution 
through "vacuum fluctuations" (viz., the net Ohu> part, where the nucleus is excited 1 htu 
and then deexcited by the same amount), a result that we will neglect for the moment 
since it doesn't affect the basic physics. Ignoring all but the monopole and dipole parts 



in eq. (28), we have 



N N 3 

PGT ^ Po(x) - ^q ■ V p (x) - j^ q ■ V q 2 V 2 p (x)\ 3hLU + ■ ■ ■ . (29) 

We perform the usual decomposition of q in terms of normalized (Cartesian) creation and 
destruction operators 

h 

(30) 



( " t+ °y5 



2\x<jj 

with 

[ai,a^] = Sij , (31) 

where /i _1 = [Zm]^ 1 + [iVm] -1 is the inverse reduced mass of the (rigid) protons-neutrons 
system, and m is the nucleon mass. 



We discuss the effects of vacuum fluctuations under ref. [17] in the reference list. The 
only modification they produce is the replacement of po in eq. (|29|) by the complete ground- 
state charge density of the model, p c ^. This "renormalization" (removal of the vacuum 
fluctuations) has for example shifted the lfiuj component of the operator q 2 q (the last term 
in eq. (p9|)) to the second term of that equation, and only true 3hui excitations remain from 
that operator. The modified eq. (^9|) then expresses the nuclear charge-density operator 
in terms of the ground-state charge density, the lhto transition charge density, the 3fno 
(dipole) transition density, etc. 

The contribution of the lhuj excitations to the nuclear moments do, and Oq are now 
easy to obtain. To evaluate the matrix element of the Schiff operator in eq. (|27| ) we need 
these two moments of the transition charge density, defined by the second term in eq. ( p9| ) 
(pdi P (x) = —[N/A] q- Vp c h(x) when vacuum fluctuations are included), expressed in terms 
of q, which contains the nuclear raising and lowering operators. One determines the lhuj 
components of these moments from eqs. (8-11): 

d GT = ™? , (32) 

O GT (lM = \d^ (r 2 ) ch . (33) 

The electron- nucleus coupling is defined by A/i(r) in eq. (|lq). Recalling that p mon in that 
equation is just p c h here, and using eq. ( |32"| ) in pd ip , we immediately obtain 

Ah GT (lhuj) = , (34) 

implying that 

5 G t(1M = . (35) 



The last conclusion also follows directly from eqs. ( p3[ ) and ([171). Equation (p4|) is more 
general than eq. (|35|), however, because it is true to all orders in Rn/Ra- 

Thus, despite the fact that it contains all the dipole strength, the Goldhaber- Teller 
giant resonance generates no contribution to the Schiff moment and therefore to the atomic 
dipole moment. In the next section we will see that the same is nearly true in real nuclei. 

4 Schiff-strength distributions, octupole correlations, and 
Schiff moments in light nuclei. 

To understand in greater detail the distribution of Schiff strength and the resulting 
ground-state Schiff moment, we first examine the situation in light nuclei. Although the 
small radii and charges of light nuclei mean that their Schiff moments will not be large 
compared to those of heavy nuclei, they have the advantage that their structure can be 
calculated at a detailed microscopic level. 

What kinds of excitations will carry the Schiff strength? We describe elementary 
excitations in terms of harmonic-oscillator shell-model quanta (hu); these are not exactly 
the same as the tuv of the Goldhaber- Teller model, but the two are related. As shown 



by Brink) 18], the electric dipole operator oIq can excite only those components of the 
harmonic-oscillator shell-model Hamiltonian (with no residual interactions) corresponding 
to the giant resonance. That is, the simple lhuj electric dipole excitations in the harmonic 
oscillator are exactly the same as in the Goldhaber- Teller model. The operator Oq, on the 
other hand, can excite other shell-model modes — isoscalar lhto and all kinds of 3hco - 
that are not a part of the GT model. The GT part of the lhuj excitations will cancel in 
S as shown above; the isoscalar lfao and all 3Huj excitations, by contrast, will contribute 
to S. 

Octupole modes are important because the E3 and Oq operators are in some sense 
the L = 3 and L = 1 angular-momentum projections of the same operator. The isoscalar 
Oq and E3 strengths, which contain both IHuj and 3hu> components, are pulled down into 
low-lying 1~ and 3~ states (in even-even nuclei) with similar structure. In combination 
with the suppression of the Schiff strength in the giant-resonance region, this similarity 
in structure results in most of the available Schiff strength being strongly correlated with 
octupole excitations. The correlation can be seen in the closed-shell nucleus 16 even 
without much calculation. The lowest 3~(6.05 MeV; T=0) state in 16 has an enhanced 
E3 transition to the ground state, B(E3)=13.5±0.7 Weisskopf units (W.u.). The lowest 
1~(7.12 MeV;T=0) state decays to this 3 _ state through an enhanced E2 transition, 
B(E2)=21 ± 5 W.u., suggesting that the two states are of similar structure (i.e. that 
the 1~ state is a quadrupole phonon coupled to the 3~ state). The 1~ state shows an 
enhanced isospin-forbidden El transition to the ground state, B(E1)=(3.6 ± 0.4) x 10 -4 
W.u. Although isospin mixing obviously contributes to the transition, a significant portion 
of the isospin-forbidden El strength appears to come from a large isoscalar matrix element 



of Oq, a part of the El operator that is normally masked [19]. 



Shell-model calculations reflect the strong correlation between the lowest isoscalar 



and 1 states [20, EIF. Furthermore, the 2p — 2h ground-state correlations in such 



calculations have a large overlap with the state formed by acting on the closed shell 
with two successive E3 operators. These higher-Tiu; correlations in the ground-state wave 
function enhance the excitation of octupole-like 3huj components in the 1~ and 3~ states, 
leading to larger low-lying isoscalar E3 and Oq matrix elements. 

Let us see how this physics works out in 19 F, which has odd A and is therefore able 
to have a ground-state Schiff moment. The ground state of 19 F is the l/2 + head of a 
K = l/2 + (sd) 3 rotational band, while the first excited state (a l/2~ state at 110 keV) is 
the basis for a 4p-lh (p _1 (sd) 4 ) K = l/2~ band. SU(3)-basis shell-model calculations p2| 
for these two bands indicate that the l/2~ band is an octupole excitation of the ground- 
state band (in fact 19 F is considered the closest thing among light nuclei|22] to an octupole- 



deformed system, though a nonnegligible part of the l/2~ state is an isovector excitation). 
Though this relation in itself implies large E3 and Schiff matrix elements between these 
two bands, we have to go beyond these old restricted calculations to get the full picture, 
for two reasons. First, the just-mentioned octupole ground-state correlations that further 
enhance the Schiff strength were omitted. Second, we need to ensure realistic behavior 
for the strength distribution of the P,T-violating NN interaction, which according to eq. 
( p7[ ) is as important as the Schiff-strength distribution in determining the ground-state 
Schiff moment. Fortunately, its behavior is much simpler. For the sake of pedagogy, we 
consider in this section a simple but often accurate[24] one-body approximation to Vpt'- 



const x a ■ f in both isovector and isoscalar channels (we will use a more sophisticated 
one-body approximation later). The lfioj El and Vpt operators then differ only in their 
effect on spin (they are in the same SU(4) multiplet). The strength from the isovector a -r 
is therefore concentrated in an SU(4) analog of the giant dipole resonance, and to first 
approximation lies at the same energy. The isoscalar strength is also mostly in a resonance 



at a similar energy, so that low- lying P,T- violating strength is depleted [25]. There will be 
some strength at low energies where the Schiff strength is concentrated, just as there is 
some El strength, but it will represent the tail of a resonance. This tail was not included 
in the calculations of ref. p^] . 

Figures 1-3 display the results of a complete (0+l)?iw shell-model calculation for 19 F, 
with the center-of-mass motion fully eliminated. In these calculations contributions to 
the Schiff moment from 2 and 2>fiu) excitations have been included through the use of an 
effective charge^]. Figure 1 shows the isovector Oq distribution, and above it the isovector 
S strength. The extent of the cancellation suggested by the Goldhaber- Teller model is 
remarkable. Also displayed is the isoscalar Schiff strength, which is uncancelled and has 
a significant low- lying component (again correlated with the E3 distribution). Figure 
2 displays the calculated a ■ r strength, and the giant resonances (including tails) are 
evident in both the isovector and isoscalar channels. Finally, in fig. 3 we graph the terms 
in eq. (^7]) as a function of excitation energy, assuming that the isoscalar and isovector 
potentials have equal strength; from this it is clear that the lowest l/2~ state almost 

4 It is difficult to treat the important 2Tiu> and ZTiui excitations consistently in the low-lying states. 
However, the concept of an effective charge for E3 transitions works well throughout this mass region. The 
results of very truncated {Q+1-\-2-\~2>)Tiuj calculations suggest a similar prescription for the Oo operator, 
and in our Yhw calculations we applied octupole effective charges to Oo- Reference fl9|] makes a case for 
this same kind of renormalization. 
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Figure 1: The distribution of isovector Oq strength (bottom), isovector Schiff strength 
(middle) and isoscalar Schiff (oc Oq) strength (top) in a shell-model calculation of 19 F. See 
text for discussion. 

completely determines the Schiff moment, which is given by the sum of all the lines in the 
plot. 

We expect the gross features of the Schiff and Vpt distributions to be general. The 
Schiff strength will be correlated with the E3 strength and lie low in energy. Nuclei with 
octupole deformation, where the E3 strength lies as low as ~50 keV and is particularly 
concentrated, will show the most-enhanced Schiff moments. Nuclei with strong low-lying 
octupole vibrations should also show enhancement. As the E3 strength moves up in energy 
(and/or the octupole collectivity is diluted), so should the dominant contribution to the 
Schiff strength, and the Schiff moments will become smaller. In the remainder of this 
paper we examine the extent to which these statements are true in heavy nuclei, where 
shell-model calculations are not possible. 
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Figure 2: Distributions of the isovector (bottom) and isoscalar (top) a ■ r strengths in 19 F. 
The corresponding operators are approximations to the P,T-violating nuclear potential. 

5 Collective Schiff moments in heavy oct up ole- deformed nu- 
clei. 

We begin our discussion with a short review of the arguments of ref. Q, the crux of 
which is related to what we have already discussed: the ground-state Schiff moment need 
not be directly related to the dipole moment, and in the case of octupole-deformed nuclei is 
considerably more enhanced. The authors adopted the particle-rotor model, which is not 
fully microscopic and omits a certain amount of valence-space physics. The arguments are 
based on collective octupole correlations, however, and this model represents them clearly 
and efficiently. 

In the particle-rotor model the nucleus is described as a single particle coupled to a 
collective core, the shape of which can be specified through a function that describes the 
dependence of its radius on angle: 



R(9, 



Bo i + B- 1 )^ 



(36) 
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Figure 3: The contributions of individual states in 19 F to the ground-state Schiff moment 
through eq. (p7|). The first excited state is dominant. 

where Rq = 1.2A 1 ' 3 . In the intrinsic frame of a deformed nucleus three of the a variables 
are no longer independent; they are replaced by three Euler angles 0{ (i = 1,2,3) that 
specify the orientation with respect to the laboratory. We will assume axial symmetry 
so that all intrinsic a's vanish except for the a/,o' s , which we denote by fy. The valence 
particle, in the "strong-coupling" version of the model, moves in a potential that is de- 
formed to match the shape of the core. The full nuclear wave function (for a state with 
angular momentum quantum numbers J and M, intrinsic magnetic quantum number K, 
and parity p) thus depends on the Euler angles, the intrinsic deformation parameters f3i, 
and the intrinsic space (r) and spin (s) coordinates of the odd particle: 



,Pl,r, s\^jmk,p) =M I + R2 



D 



J* 1 



1+pP ({3i,r,s\i& 



int/ 



(37) 



where the intrinsic particle-core state factorizes asf] 

(fl,r,s |fr int ) = $(A)V>K(r, a) . (38) 

5 In axially symmetric nuclei, the component K of total angular momentum along the symmetry axis is 
equal to that of the particle's angular momentum (sometimes denoted fi) . 
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Here M is a normalization constant, R2 rotates the /^-functions and the intrinsic wave 
function by 180 degrees around the y-axis, and P changes the sign of r and of the odd- 
multipole /3/'s in the intrinsic wave function. As long as we do not allow nonaxial core 
vibrations to be excited, we can, for the purposes of this paper, write the Schiff operator 
S z as 

S^DlKe^tCPur^) , (39) 

where S- m % is the intrinsic-frame operator, with only the K = component relevant because 
others do not generate collective excitations in the absence of nonaxial deformation or 
vibration. [For Vpt, the transformation to the intrinsic frame is trivial because that 
operator is invariant under rotation.] We have placed hats on the /3's and (on S{ n t) 
because at the quantum level they are operators that act on the wave functions $(/?;) in 
eq. (pq). All of this formalism can be justified at least in part through projected mean-field 
calculations, in which the state is a function of all A nucleon coordinates; we will argue 
shortly that such calculations are necessary to answer questions that arise in the simpler 
description. 

In the simple particle-rotor model, matrix elements of one-body operators like S z 
are straightforward to calculate. The coordinates of the individual nucleons in the core 
are integrated over in the intrinsic frame to give an operator that is a function of the 
collective coordinates and those of the odd particle, and can be applied to wave functions 
of these same coordinates. When one does the integration for the Schiff operator, assuming 
a nuclear charge density proportional to the mass density (so that the intrinsic dipole 
moment is zero^j), the result is 

§ iQt = ZeRl— Y A + ^y'+i + s int (s.p.) , (40) 

°20vr^ A /(2Z + l)(2/ + 3) 

where the last term is the contribution of the valence particle, which in this case can be 
neglected compared to the collective piece. Without the hats, the first term in eq. ( f40| ) is 
just the classical Schiff moment of a deformed drop. 

Interesting things happen when a core is both quadrupole and octupole deformed (i.e. 
when <3?(/3;) is peaked around nonzero values of 02 and (3%). [For a comprehensive review 



of the subject, see ref. [23].] The reflection asymmetry implies a double-well potential in 
the coordinate fy, which in turn means that the wave functions with good parity will be 
linear combinations of functions peaked around some (3% and its negative. If the barrier 
between the two wells is high enough, the result will be parity-doubling; low-lying states 
will have partners nearby with the same angular momentum but opposite parity. This 
means, for one thing, that there will be an intermediate state in ( BTO that enters with a 
small energy denominator and therefore a large amplitude. In fact, it is often reasonable 
to ignore all other states in the sum (see the shell- model result for 19 F in fig. 3), mainly 
because of the energy denominator, but also because if the deformation is strong enough 
the matrix element of the Schiff operator to that state is likely to be large. The reason 



6 The collective contribution to the dipole moment is obviously hard to calculate precisely if in this 
approximation it is identically zero. 
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for that is that the doublets can be viewed as the projection onto positive and negative 
parity of the same reflection- asymmetric intrinsic state |^i n t)- 

The existence of more than one state with the same intrinsic structure is exactly the 
same phenomenon as the rotational bands associated with ordinary reflection-symmetric 
deformation. The matrix element of an operator between the two states of the same dou- 
blet is proportional to the diagonal intrinsic-state matrix element of the intrinsic operator, 
just as it is for states within a rotational band. For the operator S{ n t, this diagonal matrix 
element is large; it is given roughly by the expression in eq. (^0|) with the operators 0i 
replaced by the values around which the wave function <&($) is peaked (i.e. by the classi- 
cal Schiff moment of the deformed asymmetric core) , with coherent contributions from all 
the nucleons in ity. Thus, from eq. (f27|), the physical ground-state Schiff moment S that 
determines the atomic electric dipole moment should be approximately 

s ~ -2-^sw — xe — ' (41) 

where \^o) is the opposite-parity partner of the ground state \^o), AE is the energy differ- 
ence between the two states, and S- m t is the intrinsic Schiff moment (the J-dependent factor 
is from the Euler-angle integration). The conclusion from all this is that large intrinsic 
Schiff moments and small energy denominators should make atoms with octupole-deformed 
nuclei especially sensitive tests of P,T-violation in the nucleon-nucleon interaction. 

Just how sensitive the tests will be depends on the matrix element of the interaction 
Vpt, to which we now turn. We will supply more detail in addressing this subject because 
the treatment of ref. || is incomplete and not entirely accurate. We assume that the pion 
is responsible for transmitting most of the T-violating force from one nucleon to another 
[p4| . The nucleon-nucleon interaction produced by one-pion exchange then has the general 
form|, |U 



Vpt(12) 



(<7i - a 2 ) • (r! - r~2) [Con • f 2 + C\{j\ z + t 2z ) + C 2 (3ti z t 22 
+Ci(oi + <f 2 ) • (rl - f* 2 )(T\ z - t 2z ) 



exp^-m^ri - r 2 | 
m n \fl — r^l 2 



i + ' 



vn-n \r\ - r 2 \ 




where the Cj's label isoscalar, isovector, and isotensor contributions. As long as we excite 
no particles out of the core, P, T-violating interactions between these particles sum to zero 
(see ref. H for discussion), and we need worry only about the interactions between the 
valence particle and those in the core. In what follows, we take the total mass density 
to be proportional to the charge density, and in fact from now on use the symbol p to 
represent the mass density. This assumption means that the terms with different isospin 
structure in eq. (fJ2|) enter in similar ways. Taking the range of the pion to be very short (a 
decent approximation p5[), summing eq. ([42|) over the particles in the core, and assuming 



7 As pointed out in refs. H and M, the collective enhancement of the dipole moment is much smaller 
(zero in fact if the charge distribution is proportional to the mass distribution) because the dipole moment 
is measured from the center of mass. 
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the neutron and proton densities to be equal gives the effective nucleon-core interaction 

G 



U 



PT = r] 



2mV2 



a- V/5 , 



(43) 



where (again) p is the core mass-density operator, G is the Fermi constant, and m is the 
nucleon mass (these factors are inserted to follow convention). The dimensionless param- 
eter 77 then depends on the coupling strengths Cj of the two-body interactions and the 
isospin of the nucleus. Despite its slight dependence on nuclear structure, this parameter 
is often taken as a "heuristic" fundamental quantity. The one-body approximation given 
by eq. ([£]) is slightly different from the simpler and more phenomenological one we used 
in our discussion of light nuclei. 

The form of eq. ( [43] ) makes the Schiff moment unpleasantly sensitive to the distribution 
of spin near the nuclear surface, where V/0 is largest. The results of ref. || sometimes 
differ by factors of several from those of ref. 0, primarily because of differences in the 
valence single-particle wave function, which carries all the nuclear spin s in the particle- 
rotor model. Only a significantly more sophisticated calculation (which we advocate) will 
reduce this uncertainty. We therefore do not present our own complete "particle-rotor- 
model-with-octupole-deformation" calculations in this section, but instead use that model 
in its simplest form, together with qualitative arguments, to identify a few systematic 
effects overlooked in the existing calculations. Our estimate of the size of these effects is 
obviously uncertain, but indicates what can be expected in more sophisticated calculations. 
The new physics always tends to lessen the enhancement. 

To see the what was neglected in refs. pi ||] we follow ref. [27] and expand the density 
in the deformation parameters: 



P (f) « po(r-R(6,<l>)+Ro) 

« Po (r) - R p' (r) £ m,0 + l/2R 2 Mr)C£ A^o) 



+ 



(44) 



where now po(r) is the bare (spherical) ground-state mass density of the core, a constant 
up to radius i?o m the simplest version of the liquid-drop picture. The operator a ■ V/5 in 
UpT therefore depends on the /3's, and, contrary to the statements in ref. |J, cannot be 
broken up into pseudoscalar pieces that act separately on the core and particle. In fact, 



from eq. (44) we have 



a-Vp 



= & ■ Vp 

+ RoJ2ft 



1 + 1 



21 + 



Ks-;)^'* 



I 



21 + 1 \dr 



— + ) p' [Yi-i(r] l 



i/2^EA^/Sfp(/o,ro|,o) 



v— ( ± - f\ &[Y L+ Mh - vz (I + h±± ) m.->< L 



(45) 
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where the small square brackets indicate angular-momentum coupling. Reference || con- 
siders only the first term in this expression, which changes the parity of the single particle 
and leaves the core alone. The terms with odd powers of /J3 can do the opposite, however. 
Including these other terms is important because in the strong-coupling limit of the 
particle-rotor model they tend to cancel the first term. The reason is hinted at in ref. Q, 
where it is argued that Schiff moments between close-lying states are suppressed when 
deformation is rigid and symmetric. Following ref. ||, the authors note that to the extent 
that the density is proportional to the strong one-body potential ?7 s trong felt by the odd 
particle and that the spin-orbit force is negligible, matrix elements of a ■ Vp between two 
states should be proportional to the energy difference between those states, and therefore 
very small for close-lying doublets. The reason is that under these circumstances 

a ■ Vp OC <T • W s trong = *[<? • P, ^strong] = «'[<? ■ P, H(s.p.)\ , (46) 

so that for two opposite-parity states labeled a and b with the same core structure, 

(VjMKa,\UpT\* b jMK,-p) « (VvkI*? • Vp\ip b>K ) oce a -e b , (47) 

where the e's are single-particle energies. The authors then argue that complications 
associated with asymmetric deformation eliminate this effect, but in the strong-coupling 
limit the situation is even worse because now the two states \^°jmk p) an ^ \^ b jMK -») nave 
the same intrinsic structure. As mentioned above, the matrix element of any operator 
between two such states is proportional in the strong-coupling limit to the intrinsic-state 
expectation value of the operator. In nuclei with strong octupole deformation, we therefore 
have 

(* |Vpt|*0) — ► <*int I UPT | *int) = ^TT" /f (M* " Vp|^Jf) , (48) 

which, according to the argument above, should vanish. The estimates of refs. ||, || 
apparently neglect the terms in p containing the /5's that make the shapes of the density 
distribution and the potential similar, and so do not take this effect into account. 

Of course the spin-orbit force is not negligible and the intrinsic density is not exactly 
proportional to the mean field, so the cancellation will not be complete. To get a handle 
on how much the terms containing the /3's affect the matrix element of Upr, we consider 
the ratio of the matrix element in eq. ( |4"8| ) with the terms included to that without them 
(the latter being the a simplified version of the quantity calculated in refs. [||, |(§) for a large 
number of single-particle orbits. We use a deformed harmonic oscillator as a potential, 

V(r) = -mur 2 J2PlYi,o, (49) 

l 

with deformations 02, /%, and (3^ equal to those from Table I of ref. || (we ignore higher 
multipoles and neglect pairing). We take the density to be constant inside the liquid 
drop and zero outside, a distribution that has the same angular shape as the potential, 
but a significantly different radial dependence. Figure 4 shows the absolute value of the 
ratio of {4>k\& " ^pI^k) to the same matrix element without the /3-dependent terms, for 
all K = 1/2 and K = 3/2 single-particle levels in 225 Ra below 8 Uui of single-particle 
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energy. The new terms generally have the opposite sign from that of the /3-independent 
term, and the sum is usually less than the first term by itself. The average cancellation 
is less for j = 1/2 states than for states with larger spin. Whether the decrease will be 
stronger or weaker in realistic calculations is an open question. In self-consistent mean- 
field calculations, though, there is obviously a correlation between the density and the 
spin-independent part of the field, and we could well see a significant effectfl. 

100 
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Figure 4: The ratio \a- Vp| to \a- Vpo\ (see text) for all K = 1/2 (boxes) and 3/2 (crosses) 
states in a deformed asymmetric potential for 225 Ra. The ratio is usually less than one. 

Other physics neglected in the earlier work will also have an effect. The first term in 
eq. ( fi"5| ) (<r- Vpo = {l/ r )p'o& '0 1S the spin-flip analog of the electric dipole operator (with 
a difference only in radial form and isospin). We refer to fig. 2, where the strength of the 
related operator a-f in the isovector and isoscalar channels is plotted for 19 F. As noted in 
section III, the strength is clearly concentrated in resonances at about the same energy as 



We should add that the decoupling of the particle and core[B3[, taken into account neither here nor in 
ref. |H|, could act like the spin-orbit interaction to mitigate the suppression. 
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the giant dipole resonance. A simple argument^] with a schematic residual interaction 
in RPA shows that the existence of the giant dipole resonance of the core suppresses El 
transitions between low-lying single-particle states by a factor of 3 or 4 that depends only 
on the energy of the resonance and the energy at which strength would be centered if there 
were no resonance. The low-lying transitions induced by a ■ v po should be suppressed by 
roughly the same amount because the operator a ■ r is so much like r (again, see fig. 2). 
In many-body perturbation theory, this effect can be understood as core polarization: the 
residual strong interaction can create a collective 0~ p-h pair, at the same time changing 
the parity of the valence particle. The contribution to the Schiff moment coming from the 
annihilation of the pair by Upr nearly cancels the contribution coming from the direct 
action of Upt on the valence particle itself. In any event, the resonances were completely 
neglected in refs.||, ||, and the single-particle matrix elements in those papers should 
therefore probably be three or four times smaller. 

The very large Schiff moments may be saved, however, by the combination of the two 
new effects, even though each reduces Schiff moments when added to the calculations of 
refs. J5], |6[ in isolation. The spin-flip giant resonance should not affect those parts of Upt 
that contain fa and operators like [^2c]o that do not change the parity of the valence 
particle. The suppression by the residual interaction of one part of the Schiff moment by a 
factor of three without any effect on another part may upset any balance between the two 
at the mean-field level produced by the similarity between the density and the potential. 
We will need accurate microscopic calculations to test the existence of both effects and 
the extent to which they offset one another. 

Such calculations are within the range of today's mean-field technology, and in even- 
even nuclei they have already been carried outp9|, |30|], confirming the large intrinsic Schiff 



moments in the radium isotopes[31|. In self-consistent (e.g. Skyrme-HFB) calculations the 
effects of the residual interaction on the ground state are minimized. Therefore, not only 
is the relationship between the mean field and the density likely to be most accurate in this 
case, but it is also least likely to be vitiated by corrections to the single-particle picture. 
In odd nuclei, the first-order core polarization, which has to be treated as a correction to 
particle-rotor /Nilsson models, is built into the mean- field; 1-particle-l-hole excitations of 
the core do not mix with the ground state. This fact, together with a realistic two-body 
interaction that contains all multipole-multipole terms, has implications for corrections 
to low-lying transitions from resonances, as well as for the density. The incorporation 
of core polarization means that the interplay between collective excitations and low-lying 
states is already apparent at the mean-field level, or in other words that the usual first- 
order particle-phonon mixing that reduces low-lying single-particle transitions need not 
be treated by other means (e.g. the RPA). We should therefore be able at the mean-field 
level to go a long way towards quantifying the influence of shape and resonances on Schiff 
moments in octupole-deformed nuclei. 

6 Collective Schiff moments from octupole vibrations 

The question of whether the light actinides are octupole deformed has a long history. In 
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fact the question is not entirely physical — it's really about the economy of one collective- 
model basis versus another — and it should not matter so much whether the low- lying 
states in a nucleus are best described as the rotation of an octupole-deformed shape or 
as a strong low-lying octupole vibration around a rotating quadrupole shape. Collective 
Schiff moments arise in either scheme. This fact should not be surprising in light of our 
calculations in 19 F. To see how it falls out of the collective picture, we assume that the 
nuclear core has no equilibrium octupole deformation (i.e. (0%) = 0) and write the operator 
03 in terms of creation and annihilation operators: 

03 oc 6 f + b , (50) 

where b* creates an octupole phonon with (intrinsic) magnetic quantum number K = 0. It 
is then clear from eq. (^) that the Schiff operator acting on a quadrupole-deformed state 
with no octupole phonons will create an excited state with one phonon. The terms in Upt 
(see eq. ((45|)) that are proportional to 03 can then destroy the phonon, reconnecting the 
one-phonon state to the ground state and generating a collective Schiff moment through 



eq. (]27|). 

To see how big such a moment would be we need to know the matrix element of 03 
between states with zero and one phonons. As can be seen from eq. (|50|), this quantity 



is just the zero-point root-mean-square deformation, ( J {0 2 )), which we will call 0%. [In 
other words, 03 measures the spread in 0$ of the intrinsic core wave function $(/?/).] 
This quantity can be estimated from the collective (vibrational) B(E3) transition in an 
even-even neighbor. Using eq. ([44|) to lowest order in 03 one finds f32fl 



B(E3) 0+ ^ 3 - = (3/47r) 2 (ZeR 3 o ) 2 0l . (51) 

The important point is that if a collective vibration is soft the r.m.s. deformation 03 can be 
as large as the value around which the wave- function is peaked in octupole-deformed nuclei, 
and the intrinsic Schiff moment can therefore be just as large as well. In the laboratory 
(physical) Schiff moment, there is an additional factor of 03 coming from the annihilation 
of the phonon by Upt, so that naively we expect the moment to depend on the deformation 
parameters in the combination 020 2 , where an unbarred is the value around which the 
deformed wave function is peaked. The relevant quantity for octupole-deformed nuclei 
is 02 0$ (see ref. § for a discussion of why), so that if the r.m.s. octupole deformation 
03 in a vibrational nucleus is comparable to the static value 03 of the deformation in 
an octupole-deformed nucleus, any differences in Schiff moments come from the energy 
denominator, single-particle structure, or other core excitations, not from the difference 
between deformation and vibration. We will refine this statement shortly. 

First, however, we note that the terms in Upt that don't contain 03 are usually even 
more important than those just discussed, even though they don't alter the number of 
phonons, because the zero- and one-phonon states mix through the residual strong particle- 
core interaction. The approximate form of this coupling can be derived in many ways; one 
is to examine the change in energy under a small deformation of the core. Not surprisingly, 
for an oscillator single-particle potential this leads to the same interaction that appears in 
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the octupole-deformed potential of the strong-coupling scheme (see eq. (^): 

Koupl = -mw 2 ftr 2 7 3 ,o , (52) 

where u is the oscillator energy of the (symmetric) potential, /% acts on the core, and 
r 2 Y^ t Q acts on the particle. Denoting a state with n phonons and a particle in orbit ip^R 
by \n, ^ p aK ), we have for the matrix element of the interaction V^oupi between excited 
states with one phonon and the unperturbed ground state (assuming just for illustration 
that the ground state has positive parity): 

(l,^>|Kou P i|0,^) = -mu; 2 h{^ K \r 2 Y z ^l K ) . (53) 

With a value for (3% from an appropriate B(E3), we can use the "intermediate" -coupling 
scheme of ref. [33] to diagonalize H = H(s.~p.) + H (phonon) + V conp \ separately in positive- 



and negative-parity bases (H (phonon) just contains the diagonal vibrational energies of 
the zero- and one-phonon states), so that the ground state has the form 

|* > =E^I°><K> +E^'l 1 <^>> , (54) 

8 3 

and the excited states of opposite parity have the form 

|^) = ECM|l><K)+EA,i|0,^>) , (55) 

i J 

where the i/ifx and ip7 K label single-particle states around the Fermi surface, and we are 
still ignoring nonaxial vibrations. The terms in Upt that are independent of 0$ connect 
the first terms in eq. (p4j) to the second in eq. (|55|) and vice versa. The Schiff operator 
affects the core, connecting the first term in eq. ( |54|) to the first in eq. (|55|), and the 
second to the second, effectively replacing $3 in eq. (|40|) by ^3. In this way the spherical 
/33-independent part of Upt (the only part considered in refs. || |6|) can also generate a 
collective Schiff moment. 

It is possible to use the intermediate-coupling scheme even as the phonon energy goes 
to zero and octupole deformation sets in. In that case, because the single-particle Hamil- 
tonians in the two schemes are the same, energies and matrix elements should not depend 
strongly on which scheme is used[]. One implication (which is a stronger version of a re- 
mark made above) is that if the dynamic fy associated with the vibration is comparable 
to the static /% in an octupole-deformed nucleus, and if the energy of the octupole phonon 
is small compared to typical single-particle splittings or nonaxial core-excitation energies, 
the only major difference between Schiff moments in the two cases is the energy denom- 
inator in eq. (W^j). To see this, one can imagine treating the phonon as a "decoupling" 
perturbation (along with the Coriolis interaction) in the strong-coupling scheme, as is 
done in ref. J33|. Although the diagonal matrix elements of the perturbation cause energy 



9 They will not be identical because in the intermediate-coupling scheme some of the states are particles 
and some holes, and single-particle excitation energies are measured with respect to the Fermi surface [p3|. 
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shifts, wave functions are only affected by the off-diagonal matrix elements with bands 
built on higher single-particle states or other kinds of vibrations. Thus wave functions, 
transition amplitudes, etc., will not undergo large changes until the energy of the phonon 
approaches those of other excitations. For the intermediate-coupling scheme, this means 
that the matrix elements of S and Vpt connecting ground states to low-energy octupole 
phonons should not undergo radical change once the phonon is low enough in energy, and 
that nothing special will happen in the limit that the phonons have zero energy and the 
core develops static deformation. Of course if the phonon lies high in the spectrum, the 
matrix elements can be very different from the static limit, and one must carry out the 
intermediate-coupling calculation to get a handle on the size of the Schiff moment induced 
by vibrations. 

We have done just that in several quadrupole-deformed nuclei, taking vibrational ^3's 



and phonon energies from tabulations of nearby even-even nuclei [ 32 1, and again neglecting 
pairing. In 199 Hg, the most accurately measured isotope at present, we use 03 = .09, a 
phonon energy of 3 MeV (both taken from an E3 transition in 204 Hg, which may have a 
larger Ps than 199 Hg), and quadrupole and hexadecupole parameters from ref. |jj] (for this 
simple estimate we ignore the fact that this nucleus is probably very soft). The resulting 
Schiff moment is 8 x ICP 9 ?? e f m , about half of the estimate from ref. p5[ that includes 
no nuclear correlations of any kind. In a nucleus like this, moreover, with a relatively 
high-energy phonon, the nonaxial octupole vibrations will lie nearby in energy and can 
be expected to contribute comparable amounts. When all is said and done, vibrations 
may turn out to be the dominant contribution to the Schiff moment in 199 Hg, and they 
clearly should be included in any realistic calculation. Such a calculation has never been 
done, but is crucial if we want a reliable assessment of the advantages offered by nuclei 
with strong octupole correlations. Here we need a good microscopic treatment of all kinds 
of vibrations, including the very soft 7 quadrupole mode, and must obviously go beyond 
mean-field theory. A shell-model calculation may be possible. 

An example of a large vibrational Schiff moment is in the nucleus 239 Pu. This isotope 
has several features that make it attractive for experiment, in particular a spin- 1/2 ground 
state (to eliminate quadrupole effects in a magnetic field) and a long half-life compared to 
the light actinides [The drawback is in its electronic structure, which is more complicated 
than that of Radium]. The collective E3 in 238 Pu gives /% = .09, and the nucleus has large 
quadrupole and hexadecupole deformations (J3% = .223, /?4 = .095). Together with the high 
value of Z, this makes the intrinsic Schiff moment very large. The phonon lies at 470 keV, 
about 8 times higher than the lowest state in 225 Ra, but the large intrinsic Schiff moment 
compensates in part. Our calculations, with only the /3-independent terms included in p 
(as in ref. §) give a laboratory Schiff moment of 7 x 10 7 ?7 e f m , a value a few times 
smaller than the results of ref. |§ for most of the light actinides. When we include the /in- 
dependent terms, this number goes up to 4 x 10~ 6 r/ e fro. , which is 300 times the estimate 



for Hg in ref. [35] and comparable to the results of ref. || for 225 Ra. These calculations are 
far from perfect; we had to push the energy of the octupole phonon well above the value 
from 238 Pu to get the energy of the first excited state in 239 Pu right. The uncertainty 
in the results is therefore quite large and we need microscopic calculations here too. But 
the intrinsic Schiff moments will of nuclei with low-lying octupole vibrations will clearly 
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be collective, and some may of these nuclides may be easier to investigate experimentally 
than the short-lived Radium isotopes. 

7 Conclusion 

The size of Schiff moments in nuclei with octupole correlations is determined by three 
factors: intrinsic Schiff moments, energy denominators, and the matrix elements of Vpt- 
In their discussion of the first two of these, refs. j|, g] are on rather firm ground; it is hard 
to imagine, for example, that the matrix elements between parity doublets of the Schiff 
operator are radically different from those estimates, and as we have pointed out, even 
nuclei without asymmetrically deformed cores can benefit from the same mechanism. The 
third factor is far trickier, however. 

The particle-core calculations reported both here and in earlier work can only supply 
a gross estimate of the matrix element of Vpt- The mixing that that interaction induces 
depends sensitively on the valence single-particle wave function at the nuclear surface, 
where Vp is largest. Truly microscopic calculations will give better valence wave functions 
and, if they are self-consistent, will also better represent the correlation between density 
and mean field, and incorporate the effects of resonances caused by the residual interaction. 
In vibrational nuclei it will be necessary to go a little further, but even there mean-field 
calculations will shed light on the issues we've discussed. 

Finally, for experimentalists to draw strong conclusions about enhancements over 
199 Hg, better calculations in that nucleus must be done as well. It is conceivable that 
the atomic dipole moment of 225 Ra is 400 or more times larger than that of 199 Hg (this 
is the figure reported in ref. ||), but we have pointed to physical effects that could make 
the Schiff moment in 199 Hg a few times larger than earlier calculations indicate and the 
Schiff moments in the light actinides somewhat smaller than suggested by the calculations 
of refs. [||, |(|, even within the same model. The machinery of modern nuclear structure 
theory, which is powerful enough to provide reasonably accurate estimates of the moments 
in both kinds of nuclei, should be used as soon as possible to provide experimentalists firm 
predictions for the enhancement they can expect in difficult experiments with radioactive 
nuclei. 
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